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Abstract — The paper explores the use of open-celled metal foams as compact heat exchangers, exploiting 
convective cooling. An analytical model is developed for model foams with simple cubic unit cells consist- 
ing of heated slender cylinders, based on existing heat transfer data on convective crossflow through cylin- 
der banks. A foam-filled channel having constant wall temperatures is analyzed to obtain the temperature 
distribution inside the channel as a function of foam density, cell size and other pertinent heat transfer par- 
ameters. Two characteristic length scales of importance to the problem are discussed: the minimum channel 
length required for heating the fluid to its goal temperature and the thermal entry length beyond which the 
transfer of heat between fluid and channel wall assumes a constant coefficient. The overall heat transfer 
coefficient of the heat exchanging system is calculated, and the pressure drop experienced by the fluid flow 
obtained. These results are used to analyze and guide the design of optimum foam structures that would 
maximize heat transfer per unit pumping power. Two examples are given to demonstrate the applicability 
of the analytical model: heat sinks for high power electronic devices and multi-layered heat exchangers for 
aeronautical applications. The present model perhaps oversimplifies the calculation of transport in a metal 
foam consisting of non-circular, possibly sharp-edged ligaments, and so likely leads to overestimates. 
Nevertheless the trends of heat transfer predicted by the model (for dependence on foam relative density, 
duct geometries, fluid velocity, etc.) are expected to be valid for a wide range of open-cell foams and are in 
reasonable agreement with available experimental data on aluminum foams (Bastawros and Evans, Proc. 
Symp. Application of Heat Transfer in Microelectronics Packaging. IMECE, Dallas, TX, 1997). © 1998 
Acta Melallurgica Inc. 



1. INTRODUCTION 

Non-metallic cellular foams, especially those based 
on polymers and ceramics, are widely used in ther- 
mal-insulation applications where their low thermal 
conductivity is exploited. The long list of appli- 
cation areas includes packaging of food, disposable 
hot-drink cups [1], chemical catalytic reactors, 
packed cryogenic microsphere insulations [2-4], 
solar-energy utilization, transpiration cooling [5], 
and cavity-wall insulation (dwellings, aircraft, sub- 
marine cabins, and so on). Heat is transported 
through these foams via three competing mechan- 
isms: solid conduction, thermal radiation and fluid 
conduction. (Natural) convection within the cell is 
negligible, if the cell has a size less than 10 mm [1] 
and. on the scale of a cell, the typical temperature 
differences between the fluid and the solid cell edges 
are small. At room temperature, the thermal con- 
ductivities of most non-metal foams are well below 
1 W/(m K), with those of polystyrene, polyurethane 
and phenolic foams having low relative densities 
(-0.02) roughly equal that of air, -0.025 W/(m K). 

Metal foams, on the other hand, are more likely 
to find applications in ultralight structures where 
stiffness, strength and toughness are emphasized. 
The effective thermal conductivities of open-celled 
metal foams, about one to two thirds of the pro- 
duct of their relative density and the thermal con- 



ductivity of the solid of which they are made, are at 
least one order of magnitude larger than those of 
their non-metallic counterparts [1]. hence are gener- 
ally unsuited for thermal-insulation purposes. 
Open-celled metal foams, however, can be used to 
enhance heat transfer in applications such as cryo- 
genic heat exchangers, heat exchangers for airborne 
equipment, coal combustors [2], compact heat sinks 
for high power electronic devices [6], heat shielding 
for aircraft exhaust, compact heat exchangers [7], 
liquid heat exchangers [7,8], air-cooled condenser- 
cooling towers and regenerators for thermal 
engines [4]. More uses of these relatively new ma- 
terials are expected in the future, particularly 
because low density foams with remarkably uniform 
and regular cell morphologies are currently being 
developed using affordable processing methods. 

Most of the previous studies on heat transfer in 
porous medium are based on Darcy's law which is 
only valid for creeping flows through an infinitely 
extended uniform medium with Reynolds numbers 
less than about 10 [2,9]. Under such conditions, 
the fluid and the solid matrix may be assumed to 
be in local thermal equilibrium so that the fluid 
saturated porous medium can be treated as a 
continuum [7,8,10,11]. This paper deals with 
direct-contact type heat exchangers (e.g. a porous 
channel bounded by isothermal or isoflux plates) 
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where the fluid velocity is high and the porous med- 
ium is bounded, rendering the assumption of 
Darcian flow and local thermodynamic equilibrium 
invalid. The objective is to obtain a functional re- 
lationship between the cellular structure and the 
heat-transfer characteristics for forced convective 
flows through open-celled metal foams. A simple 
cubic cell model consisting of slender cylinders as 
edges is developed to capture the most important 
behavioral trends of energy flow due to forced con- 
vection and conduction through cell ligaments of 
the cellular foam. Use is made of the analogy 
between flow through the foam and that across a 
bank of cylinders. The solid matrix and the convec- 
tive fluid are treated separately by the model. 
Approximate closed-form solutions for the overall 
heat transfer coefficient and pressure drop as func- 
tions of cell morphologies and other relevant heat 
transfer parameters are obtained, which can be used 
to guide the design of optimum cellular structures 
that maximize the heat transfer rate per unit pump- 
ing power (the energy expended driving the convec- 
tive flow). The model is applied to evaluate the heat 
transfer efficiency of open-celled foams as heat 
sinks for power electronic components, as well as 
the performance of air-fuel heat exchangers consist- 
ing of foam laminates. 

It should be pointed out that the present model, 
which is based on crossflow across a bank of cylin- 
ders, has been idealized in several ways to limit the 
complexity of forced convective flow across an 
open-celled foam. Some of the simplifying assump- 
tions made in order to analyze heat transfer in the 
disordered porous material likely lead to an over- 
prediction of the actual level of heat transfer but 
nevertheless should capture the approximate func- 
tional dependence of the different control variables 
such as foam density, cell size and fluid velocity [2], 




Fig. 1. Open-ceiled foam under forced convection: (a) 
notations, (b) cubic unit cell. 



2. SPECIFICATION OF THE MODEL 

Referring to Fig. 1(a), we consider the steady- 
state heat flow in an open-celled foam made up of 
uniform distributed, equal-sized cubic cells, sand- 
wiched between two plates of length L and width 
W. The cubic unit cell is chosen for its simplicity, 
allowing for approximate closed-form solutions for 
important heat transfer parameters. The dimensions 
of the cell are chosen such that the surface area 
density of the model foam equals that of the real 
foam which can be measured accurately by using a 
micro-computed-tomography system [2]. Both plates 
are assumed to be thin and have large thermal con- 
ductivity so that the through-thickness heat conduc- 
tion may be neglected; thus, uniform temperatures 
T\ and T 2 are prescribed, respectively, on the top 
(z = 0) and bottom (z = b) plates. Without loss of 
generality, it is assumed that T X >T 2 , and that the 
sandwich structure is capped and thermally insu- 
lated at both ends y = Wjl and y - - Wjl. 
Cooling fluid, velocity v 0 , temperature T 0 ( < T 2 ) 
and pressure p 0 , is forced into the foam at x = 0 
(the inlet) and exits at x - L (the outlet) with tem- 
perature T e and pressure p c . The simple cubic unit 
cell shown in Fig. 1(b) is taken as consisting of 
three mutually perpendicular equivalent cylinders, 
each having diameter d, length a and thermal con- 
ductivity A s . For simplicity, the cylinders making up 
the cell edges are assumed to be each parallel to the 
x-, y- or z-coordinate and referred to hereafter as 
cylinders X, Y and Z, respectively. Thus, the flow of 
fluid is normal to cylinders of both types Y and Z, 
but parallel to type X (Fig. 1). The cell size a is 
assumed to be much smaller than the plate width W 
so that heat transfer in each cell envisioned in Fig. 1 
is assumed to be independent of y. The design objec- 
tive is to obtain the maximum heat transfer per unit 
pumping power from the sandwich system. 

The model presented in Fig. 1 is analyzed in 
three steps. We start in Section 3 by calculating the 
heat transfer from a single cylinder of type Z and 
length b y with temperature T\ specified on its top 
and T 2 at its bottom (Fig. 2(a)). When T 2 is suffi- 
ciently low, heat is lost not only by convection to 
the fluid but also by conduction to the environment 
(through the lower plate). However, for higher T 2 
or a larger heat transfer coefficient on the cylindri- 
cal surface, reverse flow of heat occurs at the bot- 
tom half of the cylinder, so that heat is extracted 
from both the lower plate and the upper one into 
the cylinder. The condition governing such reverse 
flow of heat is obtained. These results are then 
extended in Section 4 to analyze a bank of parallel, 
unconnected cylinders of type Z in cross-flow, 
arranged in a square array in the (x, y)- plane. 
Expressions for the average heat transfer coefficient, 
average fluid temperature inside the array, and total 
heat transfer from the system are obtained. The 
problem specified in Fig. I is solved in Section 5, 
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(a) 




(b) 

Fig. 2. Single cylinder in cross-flow: (a) longitudinal view, 
(b) cross-sectional view. 

based on solutions described above. The surface 
area density of the model foam is calculated and 
compared to the data provided by ERG (Energy 
Research and Generation) for their Duocel® alumi- 
num foam, as well as data for other types of heat 
exchangers. The overall heat transfer coefficient and 
pressure drop caused as the fluid flows across the 
foam are calculated using the model and the opti- 
mum foam properties for the best heat transfer per- 
formance are suggested. Finally, two examples of 
practical importance are presented in Section 6: 
compact multi-chip modules having convectively 
cooled foams as heat sinks, and multi-layered heat 
exchangers where heat transferred from one flowing 
fluid to another is substantially enhanced by the use 
of foams. 



3. A SINGLE CYLINDER IN CROSS-FLOW 

The metal ligaments in an open-celled foam are 
similar in many ways to fins common in heat 
exchangers so we begin by considering a single, iso- 
lated cylinder of type Z. originally situated at (jc, y) 
in the sandwich structure (Fig. 2). The average tem- 
perature of the fluid flowing across the cylinder. 7>, 
depends on the position of the cylinder in the direc- 
tion of flow (x-direction), but is here assumed to be 
independent of z. (The validity of this assumption 
and its limitation is discussed in Section 5.2. 
whereas the x-dependence of 7> is treated in Section 
4.2.) Let h be the average heat transfer coefficient 



associated with the cylinder subjected to cross flow. 
The cylinder is considered to be slender 
(d « a < b) so that its temperature, T(z), is uni- 
form across the cylinder which allows for a simple 
closed-form solution for 7\z). This simplification 
stems from the observation that the aspect ratio dja 
is typically small for most low-density, open-celled 
metal foams (Duocel & Al foams, for i nstance, have 
an aspect ratio ranging from 0.1 to 0.2, see Fig. 4). 

S.I. Temperature field and heat transfer 

A standard "fin" analysis is useful for character- 
ization of the metal ligaments. Under steady-state 
conditions, the variation of temperature T along the 
length of a slender cylinder is governed by [13] 



d 2 T Ah , 

___ (r _ rf) = 0 



(3-D 



where the Fourier law of heat conduction and 
Newton's law of cooling have been used. For typi- 
cal uses of metal foams, heat transfer due to radi- 
ation at the cylindrical surface is usually at least 
one order of magnitude smaller than that due to 
forced convection (Appendix A), so it will be neg- 
lected. The solution to equation (3.1), with the 
boundary conditions that at r= T 2 z = 0 and at 
T[ = T 2 z - b, is 



r(z) = r f + 



(7, - r f )sinh[2Bi 1/2 (A - z)jd\ 
+(r 2 ~r f )sinh(2Bi 1/2 z/^) 
sinh(2Bi ,/2 z/d) 



(3.2) 



where Bi - hd/* s is the Biot number. The average 
temperature of the cylinder, f, is therefore 



T= 7} + (r, +T 2 -27»- 



cosh(2Bi l/2 6/</)- 1 
(2Bi 1/2 £/J)sinh(2£/'/ 2 6A/) 



(3.3) 

which has a lower bound T= 7> when Bi — oo 
(convection dominated) and an upper bound 
T={T X + T 2 )I2 at Bi = 0 (insulated cylinder). In 
most applications involving air and water cooling 
A- 10 2 -l0 4 W/(m 2 K), rf-10^m and ^10 2 W/ 
(m K) rendering Bi « 1, such that (7*, + T 2 )/2 is a 
good estimate for T. Equation (3.2) dictates that a 
thermal boundary layer with thickness 0(Bi if2 d) 
exists near both ends of the cylinder for the high Bi 
regime. The temperature is more uniform away 
from the ends, approaching 7} for Bi > 0.02. In gen- 
eral. T*(T X + T 2 )I2 if Bi < 0.002 whilst f * 7> 
when Bi>0.02. 

The heat flux entering the cylinder at 
r = 0. q = -(7t</ 2 /4)A s (dr/dz) r -o. Allows from 
equation (3.2) as 

?=(^/2)>l s Bi 1 ^ 

(7j - 7Y)cosh(2Bi 1/2 fr/<y) - (7% - 7V) 
sinh(2Bi l/2 ^) 



(3.4a) 
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whereas the loss of heat at the end of the cylinder 
due to conduction through the solid, q c = -(nd 2 j 
4)>i s (dr/dz) r » b , is given by 

q c =(7t/2)A s Bi ,/2 </ • 

(ri-r f )-(r 2 -r f )cosh(2Bi 1/2 ^) 

X sinh(2Bi 1/2 *A/) 

The heat carried away by the cooling fluid, q h , is 
obtained by subtracting q c from q, yielding 

q h = (7i/2)/ s Bi 1/2 </(r, + T 2 - 2r r )tanh(2Bi ,/2 ^) 

(3.4c) 

However, depending on the magnitude of the Biot 
number Bi, q c of equation (3.4b) may be negative, 
as exploited below, owing to the competition 
between conductive loss at z = b (T\ > T 2 ) and 
positive conductive flux (T 2 > 7». 

3.2. Condition for reverse flow of heat 

In the limit when the cylinder is insulated, 
equations (3.4aHc) reduce to 

q = q c = (nd 2 l*)X % {T { - T 2 )/b, <? h = 0 (3.5) 

In the limit Bi — > oo when convective transport is 
dominant, 

q c /q = -(T 2 -T ( )/(Ti-Tt), 

q h /q = (T { + T 2 -2Tt)/(T l -Tt\ (3.6) 

so that, with T C <T 2 <T U reverse flow of heat 
occurs at the lower half of the cylinder, q c = - q 
and <7h = 2(? if 7j = T 2 . For finite Bi, heat flux into 
the system from the lower boundary at z = b 
occurs when Bi > Bi* where 

It is clear that energy flows into the cylinder from 
both ends if Tf<T 2 =T u i.e. Bi* - 0. Thus, the 
total transfer of energy to the fluid from the single 
cylinder system, q, is 

_ \q ifBi<Bi* nR v 
*=U if Bi^Bi* (3 8) 

As an illustration, given (T\ - T r )l(T 2 - 7» - 2, 
A s = 200 W/(mK), d = 0.5 mm and b = 1004 the 
heat transfer coefficient h must be larger than 10 W/ 
(m 2 K) if reverse flow of heat is to occur at the bot- 
tom of the cylinder. 



4. A BANK OF TYPE Z CYLINDERS IN CROSS 
FLOW 

In this section, the results obtained above for a 
single cylinder of type Z are applied to study a sim- 
plified version of the sandwich structure shown in 
Fig. 1(a) where the fluid passes a bank of type Z 
cylinders, squarely arranged in the (x, v) plane but 



unconnected by cylinders of either type X or Y. The 
effects of these connecting struts constitute the sub- 
ject treated in Section 5. 

4.1. Heat transfer coefficient h 

The transfer of energy due to high- Reynolds- 
number streaming flow past a heated single cylinder 
or a bank of heated cylinders has been a subject 
area of enormous studies, but results are largely 
limited to empirical correlations. For most practical 
applications, dimensional analysis of the heat trans- 
fer processes suggests that the Nusselt number, 
Nu = hdlkfi is a function of the Reynolds number, 
Re = vd/Vf, and the Prandtl number, Pr - Cp/ir/Af. 
Here, v is the average velocity of the fluid far from 
the cylinder, and A ft v f , u f and c p stand for the ther- 
mal conductivity, kinematic viscosity, shear vis- 
cosity and specific heat (at constant pressure) of the 
fluid, respectively, all evaluated at the film tempera- 
ture (the arithmetic mean between the surface and 
free-stream temperatures). The Nusselt number Nu 
may be interpreted as the ratio of the actual heat 
transfer coefficient h to the conductive heat transfer 
coefficient X ( /d of the fluid. The Prandtl number Pr 
is the ratio of the kinematic viscosity to the thermal 
diffusivity, measuring the relative efficiency of the 
fluid as a conductor of momentum and energy. The 
Reynolds number Re measures the ratio of inertial 
forces to viscous forces in the flow. Consequently, 
the overall heat transfer coefficient is expected to 
increase as the Reynolds number or Prandtl num- 
ber, or both, increase. 

According to Gosse [14], for a single cylinder in 
cross flow (Fig. 2(b)) and Re > 40, the local Nusselt 
number at the front face (0 < <f> < 4tt/9) can be 
expressed as 

N«>) = h(4>)d/M = 1.05(1 - <£ 3 )Pr 1/3 Re ,/2 

Heat transfer due to convection at the back face 
(4tt/9 < <j> < n) is relatively small for small Re num- 
bers (Re < 40), but may exceeds that at the front 
face when Re becomes large. The average Nusselt 
number for the whole cylinder may then be 
described in the form 

Nu = CPr 1/3 Re\ (4.1) 

where C and n are dimensionless constants which 
has been found to be appropriate for most ordinary 
gases and liquids, except for liquid metals [9, 15, 16]. 
This expression is also useful for non-cylindrical 
shapes of objects, subjected to changes in (C, n) 
and the replacing of d by the hydraulic diameter 

For a bank of parallel but unconnected cylinders 
subjected to cross-flow, heat transfer is still gov- 
erned by the Nusselt number expressed in the form 
of equation (4.1) such that the overall heat transfer 
coefficient averaged over all the cylindrical surfaces 
may be described by [16, 17] 
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h = MC?r^(v max /v { ) n d n ~ l (4.2) 

Here, v max is the maximum fluid velocity occurring 
in the bank, and (C, n) are dependent on the 
Reynolds number and the geometric characteristics 
(e.g. relative spacing of the cylinders, numbers of 
rows of cylinders, and the flow direction). Equation 
(4.2) has been correlated with extensive experimen- 
tal data for both gases and liquids and for a num- 
ber of in-plane arrangements of cylinders [16,18- 
20]. For cylinders arranged in a square array with 
at least 10 rows in each direction such that a repea- 
table pattern of flow is established inside the bank, 
C^0.3 and n ^ 0.6 when a/d>3, rendering 
hccct 04 [16.17]. For air at atmospheric pressure 
and room temperature with v max =] m /s and for a 
bank of cylinders having diameter d-= 0.5 mm. one 
calculates from equation (4.2) that /i^150W/ 
(nr ft). If air is replaced by water with other geo- 
metric and flow conditions unchanged 
h*4x i0 4 W/(m 2 K). 

Using the properties tabulated in Holman [16] for 
common fluids (e.g. air and saturated water), we 
find that the heat transfer coefficient h depend 
weakly on temperature. Therefore, the temperature 
dependence of h is neglected in the analysis below. 

4.2. Fluid temperature 

The temperature of the cooling fluid increases 
continuously as it flows along the bank of cylinders 
and so the driving force for heat transfer is continu- 
ally diminished. We therefore wish to establish the 
fluid temperature as a function of downstream pos- 
ition .v, the plate temperatures (T u 7\) and the 
properties of the foam. 

It is not attempted below to find the detailed 
temperature field of the fluid everywhere in the 
bank, given the turbulent flow and limited infor- 
mation on the heat transfer coefficient. Instead, 
several simplifications are made to arrive at an ap- 
proximate solution for the steady-state temperature 
distribution of the fluid inside the bank. For simpli- 
city, the lower plate temperature T 2 is set equal to 
T, although this restraint can be relaxed if necess- 
ary. We assume that heat transfer between the fluid 
and the plates is governed by the same coefficient h 
for cross-flow over cylinders, namely equation (4.2). 
We further assume that the turbulent flow inside 
the foam is well mixed and so develop a one- 
dimensional model for the fluid temperature, 7}, 
with the variations of 7> with z (and y) neglected.' 
This is partially supported by the result of 
Kaviany [9] on the fully-developed flow through a 
porous channel bounded by isothermal parallel 
plates (Fig. 1): 

K_;) _ 1 - e' 2y - (1 - e - y )[e*'*-'> + e~^ b ] 
v o I -e-2>'-2(l -e-7) 2 r i 

where y = {bh/fQ l/2 is the porous medium shape 



parameter, K being the permeability of the porous 
medium and c the porosity. There is a strong depen- 
dence of the fully-developed velocity profile v(r) on 
y; as y increases, the central region containing a uni- 
form velocity distribution spreads further toward 
the plates, with the velocity variation confined to a 
very thin layer adjacent to the plates at large y. 
Also, by dimensional analysis. Vafai and Tien [11] 
showed that the boundary-layer thickness is of the 
order of y~'. Notice that for a duct filled with metal 
foam, the typical values of y range from 10 l to 10 3 . 

Consider an arbitrary slice of the sandwich struc- 
ture located at x with length Ax (Fig. 1(a)). With 
W» = Ax/a 2 denoting the total number of cylinders 
inside the slice per unit thickness, the loss of heat 
from the cylinders to the streaming fluid is /V s times 
that from a single cylinder, q h , given by 
equation (3.4c). An energy balance on the sliced 
structure shown in Fig. 1(a) gives 

mc p [T r (x + Ax) - T t (x)] = /Vrf h + g w (4.3) 

where m = p^fi is the mass flow rate per unit width 
at the entrance to the bank and p r is the fluid den- 
sity at the entrance. Also, q w is the heat flux into 
the fluid from both plates over length A.v and unit 
width, 

?w = (2^r)A[r, - T ( (x)] (4.4) 

where 

7=1- nd 2 /4a 2 (4.5) 

is the net surface area of the plate per unit length 
per unit width excluding the cross-section areas of 
the N s cylinders. Combining equation (4.3) with 
equations (4.4) and (3.4c) (with r, = T 2 ) results in 
T r (x + Ax) - Tt(x) 2 n h 

■[■^(D'-It 5 )] <«> 

It follows immediately that the characteristic length 
scale / for fluid temperature variations in the direc- 
tion of flow is given by 

(4.7) 

In the limit Ax — 0, equation (4.6) provides an 
ordinary differential equation for Tfa) which can 
be readily solved to obtain 

T ( (x) = r, - (F, - 7i)exp(-*//) (4.8) 

The average of T t {x) over length L is 

(4.9) 
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Bi=0.05; b/d =200 ; L= t T 0 /T l= 0.2 




- 3 Predicted foam temperature distribution for (a) Bi l/2 = 0.005, (b) Bi 1/2 = 0.05 (with bfd = 200, 
' L = /and r 0 /T,= 0.2). 



The significance of equations (4.7), (4.8) and (4.9) 
will be discussed below in some detail regarding 
model foams with simple cubic cells. The tempera- 
ture of each cylinder in the bank is obtained by 
substituting equation (4.8) into equation (3.2). As 
an illustration, Fig. 3 presents the distribution of 
the normalized cylinder temperature, T(x, z)jT u 
for Bi 1/2 - 0.005 and 0.05, respectively, with 
tyd = 200, TqIT x =0.2 and the foam length L 
chosen equal to the characteristic length scale /. 
Notice the presence of a boundary layer along the 
cylinders near the plate wall, which thickens as Bi 
decreases or x/L increases. 



5. FOAMS WITH OPEN CUBIC CELLS 

The analysis of heat transfer in the model foam 
structure thus far has ignored the contributions to 
heat transfer from type X and Y cylinders that con- 
nect cylinders of type Z. As illustrated schematically 
in Fig. 1(b), the three types of cylinders (which are 
assumed to have a high conductivity relative to that 
of the fluid), mutually perpendicular, are each 
arranged in a square array in their respective cross- 
planes and meet at three-way joints uniformly dis- 
tributed inside the foam. Under steady-state flow 
conditions, the average temperature of type X and 
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Fig. 4. Predicted cylinder diameter as a function of foam density compared with that measured for 
aluminum foam, for three cell sizes: a = 0.5, I and 2 mm. 



Y cylinders inside a cell located at .v is expected to 
be approximately the same as that of the type Z 
cylinder meeting them at the joint of the cell. Thus, 
all the cylinders in a single slice of the structure 
shown in Fig. 1(a) may be assumed to have an 
average temperature T(x\ coupled to the average 
fluid temperature Tfa) through equation (3.3). h 
follows that type Y cylinders, subjected to cross 
flow, have the same average heat transfer coefficient 
h as that given by equation (4.2) for type Z cylin- 
ders, also under cross flow. On the other hand, 
although the average heat transfer coefficient for 
type X cylinders under longitudinal flow is still 
described in a form similar to that given in 
equation (4.2), extra terms need to be introduced to 
account for the dependence on the cylinder 
length [14]. Such dependence is nonetheless assumed 
to be fairly weak when x> 10a (see discussion 
below) where the flow pattern is well established. 
For simplicity and also due to the lack of more 
detailed data on longitudinal flow over a bank of 
cylinders, it is assumed below that type X cylinders 
are characterized by the average heat transfer coeffi- 
cient of equation (4.2). Under these assumptions, 
the solution obtained in Section 4 for a bank of 
type Z cylinders applies as well to the model foam 
structure of Fig. 1(a), with due account of the extra 
number provided by type X and Y cylinders plus 
their contributions to the total surface area for heat 
transfer. 



5.1. Cylinder diameter and surface area density 

For an open cubic cell consisting of slender cylin- 
ders, the diameter of the cylinder is approximately 
related to the relative foam density, p, by 



«*(v/3^)' 



(5.1) 



where a is the ceil size and p = p*/p s , with p* 
denoting the foam density and p s the density of the 
solid. The resulting total surface area A per unit 
width (in the >'-direction) of a segment of foam of 
height b and length L, is approximately 



(5.2) 



The surface area density, z A s AfV with V denoting 
the volume of the foam, is an important parameter 
for the heat exchange systems: 

The predicted cylinder diameter, from 
equation (5.1) is plotted on Fig. 4 as a function of 
relative foam density, p, for three values of cell size, 
a - 0.5, 1.0 and 2.0 mm. The ERG data for 
Duocel^ Al foam are also included, corresponding 
to cell sizes of 10, 20 and 40 ppi (pores per inch), 
respectively. The mean cell size for these foams. 
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Fig 5 Predicted surface area density as a function of foam density compared with that measured for 
aluminum foam, for three cell sizes: a = 0.5, 1 and 2 mm. 



according to the ERG catalogue, are 2.0, 1.0 and 
0.5 mm, respectively. The simple cubic cell model is 
seen to be fairly representative, bearing in mind 
that for a given ppi, the nominal cell size of a real 
foam is not constant but dependent upon p, and 
also that the cell ligaments of the ERG foams are 
best described as a rod with rectangular or triangu- 
lar cross-section rather than a circular cylinder. 

In Fig. 5, the surface area density, x A , calculated 
from equation (5.3) is plotted as a function of p 
with a = 0.5, 1.0 and 2.0 mm, respectively. Again, 
the model correlates well with the ERG data. In the 
40 ppi case, however, a 25% difference between the 
prediction and measurement is present which may 
be attributed to the actual mean cell size being 
different from the nominal mean cell size used in 
the calculation, a = 0.5 mm. (Microstructural 
characterizations suggest that the Duocel® foam 
labelled by ERG as 40 ppi is actually described 
more closely by 30 ppi.) With <x A > 1000 m 2 /m 3 even 
at very low densities (p * 0.02), the conclusion 
drawn from these results is that open-celled metal 
foams well qualify as compact heat exchangers. 
Compact heat exchanger generally require 
^>700m 2 /m 3 [9,21] and are essential in appli- 
cations where the size and weight of the heat 
exchanger is constrained due to design consider- 
ations (e.g. automobiles, airborne equipments and 
air conditioners). For comparison, automotive 
radiators have * A * 600-1000 m 2 /m 3 , cryogenic heat 
exchangers 1000-2000 m 2 /m\ gas turbine rotary 



regenerators 3000-6000 m 2 /m\ and human lungs 
10 4 m 2 /m 3 [21]. 

5.2. Characteristic length scales 

5.2.1. Minimum flow length /,. As shown in Fig. 6, 
two length scales are of importance to the present 
design problem: the length /, the fluid must travel 
to raise its temperature from T G to the goal tem- 
perature T G (which is fixed by the plate temperature 
T,), and the distance l 2 from the inlet at which the 
bulk of the fluid attains a fully mixed state with 
uniform temperature and velocity except in a region 
of thickness 6 adjacent to the plate. Practically, l x 
may serve as a guideline for selecting the total 
length L of the sandwich structure to achieve the 
targeted fluid temperature at the duct exit, whereas 
l 2 dictates the range within which the current model 
is valid. 

Consider first the characteristic length scale /, for 
fluid flow. From its definition above and 
equation (4.8), it is related to the length scale / for 
fluid temperature changes (equation (4.7)) via 

Since each cubic cell effectively contains one cylin- 
der of each type (assuming a/</>3), the total num- 
ber of cylinders, N» inside a slice of length A* and 
unit width, now becomes 



N 



(5.5) 
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Fig. 6. Characteristic length scales of heat transfer in a duct filled with open-celled foam- (a) tempera- 
ture profile, (b) velocity profile. 



Since the second term in the bracket on the right 
side of equation (4.7) has a typical value of at least 
10, / can be rewritten as 



(5.6) 



37u,Bi !/2 tanh(2Bi 1/2 6/</) 

For consistency, one needs to check if at x = /, the 
average cylinder temperature f(/,) still exceeds T G . 
From equations (3.3) and (4.8), the average tem- 
perature of the cylinders as a function of x is 
obtained as 

r ' V TJl (Bi I/2 />/t/)sinh(2Bi 1/2 V^)J 

(5.7) 



x exp(-x//) 
It follows that 
TU\-t _l cosh(2Bi l/2 £/</)- l 

r(/,) ~ rc V^^^ 

(5.8) 

which is larger than T G as long as 7* G < 7Y 

Selected calculations for air reveal that, for rela- 
tively large foam densities (p>0.2), the normalized 
minimum flow length Irfa increases sharply as the 
normalized foam thickness bja increases, becoming 
linear in bja when bja is larger than about 10. On 
the other hand, l x \a is nearly independent of p 
when the foam is reasonably thick (b/a> 10). In ad- 
dition, it has been verified that IJa decreases as the 
Biot number Bi is increased. 



5.2.2. Thermal entrance length l 2 . By definition, 
the thermal entrance length, / 2 , is the longitudinal' 
length from the inlet to the location where the ther- 
mal boundary layer becomes as thick as the duct 
itself. In contrast to the thermal entrance region 
(x < l 2 ) where the shape of the transversal tempera- 
ture profile changes from one location to the next, 
in the thermally fully developed region (x>l 2 ), the 
shape (although not the amplitude) of the tempera- 
ture profile is preserved. It is expected that the ther- 
mal entrance length is closely related to the 
hydraulic entrance length l{ (i.e. the longitudinal 
distance at which the center-line velocity is 99% of 
its value for a fully developed field). Compared 
with the other characteristic length scale /,, the 
length U is more difficult to calculate owing to the 
complex and unsteady nature of fluid flow in the 
entry region. A rough estimate of / 2 is attempted 
below. 

Consider the entry region to the duct as shown in 
Fig. 6. If the duct is not filled with foams, it has 
been established [22] that the thermal entrance 
length under laminar flow conditions is approxi- 
mately given by 

/2/£*0.05Re*Pr (5.9a) 

which is of the same order of magnitude as the 
hydraulic entrance length l{ for fluids with Pr * 1 
(e.g. air): If l{ is negligible relative to 6, then A is 
also negligible [10]. When the Reynolds number 
Re6 = Vf£/v f has a value larger than about 
2000 [15,16], turbulent flow prevails in the duct, 
resulting in [16] 
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h/b = 50 



(5.9b) 



For a duct filled with open-celled foam subjected to 
convective flows with Re 6 >2000, the fluid flow is 
expected to be more turbulent and unstable once it 
enters the duct because of the large disturbances 
caused by the presence of numerous foam cells and 
their torturous ligaments. In addition, the rough- 
ness of the plate surfaces due to the foam attach- 
ments is expected to dominate the turbulent motion 
in the boundary layer adjacent to the duct wall. 
Hence, a smaller h\b than that given by 
equations (5.9a)-(b) is expected throughout the 
duct. Kaviany [8] found that the hydraulic entrance 
length l{ decreases rapidly with an increase in the 
porous medium shape parameter -/-(be/K) 
(/V oc y~ 2 for flow over a plate in saturated porous 
media [11]). Thus, as 10 1 <y < 10 3 for a typical 
duct filled with metal foam, the hydraulic entrance 
length \{ for the porous medium is much shorter 
than that for a nonporous medium and is generally 
negligible [8]. Similar trends are expected to hold 
for the thermal entrance length / 2 . 

5.3. Overall heat transfer coefficient n 

If the amount of heat transferred directly from 
the bounding plates of the duct to the fluid, <? w , is 
neglected based on the same argument as that lead- 
ing to equation (5.6), the loss of heat from a single 
slice of foam, length A.v and unit width, would be 
N s times that from a single cylinder: g s = N s q h . The 
combination of equations (3.4c), (5.1) and (5.5) 
then leads to 

Q 5 (x) =[273^ s Bi l/ y /2 tanh^Bi 1 ^)] 

x [r, - Tr(x)]{Ax/a) (5.10) 

The total heat transfer from the whole sandwich 
structure per unit width, g s (x), is obtained by inte- 
grating over all the slices, yielding 

Q 



such that heat transfer from the system is equival- 
ent to that carried away by a fluid of average tem- 
perature f f , which flows across a foam-free but 
otherwise identical duct with the same average heat 
transfer coefficient Ji at the wall. The factor 2 is 
needed in equation (5.14a) because the duct has two 
plates each of length L. A straightforward manipu- 
lation using equations (4.8) and (4.9) shows that the 
temperature difference, T\ - f f , is the well-known 
logarithmic mean temperature difference [16], A7 m , 
namely 

T e -T Q 



A S (T, - To) 



Ar m = r, - f f = 



in[(Ti - r 0 )/(r t - r e )l 



(5.14b) 



=(l/a)[\ -exp(-L//)] 

x[2V3^Bi 1/2 p 1/2 tanh(2Bi l/2 6/^)J5.11) 

Substitution of equations (4.8), (5.1) and (5.6) into 
equation (5.11) results in 

Q = m c,(T c -T 0 ) (5.12) 

where T c is the temperature of the fluid at the outlet 
x = L (Fig. 1(a)), given by 

r c = r, -(r, - r 0 )ex P (-L//) (5.13) 

Obviously, the maximum amount of heat that 
the fluid can extract from the system is 

Qm™ = ™p(T\ ~ r o)- 

The overall heat transfer coefficient of a foam- 
filled duct, k may be defined as 



From equations (4.9), (5.11) and (5.14aHb), one 
has 

^ = (2^Br l/2 )tanh(2Bi ,/2 ^/c/) (5.15) 
which, in the limit Bi 1/2 6/</ « 1, reduces to 

h = (x A b/2)h, B\ l/2 b/d«\ (5.16a) 

In the other limit when Bi 1/2 6/</>2 such that 
tanh(2Bi I/2 6/J) * 1, 

h = (* A mJhkTd, Bi 1/2 fc/</>2 (5.16b) 

where i A is the surface area density of the foam 
given in equation (5.3). Notice that (ct A b/2)2L = A 
is the total surface area of heat transfer in the foam 
per unit width. Equation (5.15) is in general agree- 
ment with the test data of Bastawros and Evans [1] 
on ERG foams, giving the correct functional depen- 
dence of h on the fundamental variables (foam den- 
sity, channel width, fluid velocity, etc.). However, 
the predicted value of h from equation (5.15) is 
much larger than that measured: At high air vel- 
ocities (v f > I m/s), the model predicts a typical 
10-20-fold increase in h if an empty duct is filled 
with an ERG foam; a 3 ~ five-fold increase of Ji is 
measured in stead. The overestimates may be 
attributed to the simplifying assumptions of the 
model in order to obtain simple analytical solutions 
to the difficult, relatively unstudied problem. On the 
other hand, the test set-up of Bastawros and 
Evans [2] is not quite the same as that shown in 
Fig. I and it is not clear if the forced air flow has a 
uniform velocity before and after it enters the foam 
(which is assumed to be the case in the model). 
Obviously, more measurements are needed to 
further check the validity and accuracy of the 
model predictions. 

In analogy to standard duct flow problems, an 
average Nusselt number for the foam-filled duct 
system, Nu, may be defined by [16] 

m = hD/M (5.17a) 



Q = 2Lh(T x - f f ) (5.14a) ^ . g ±q hydraulic diameter, 
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n _ A x Cross-sectional area of fluid flow local Reynolds number Re in equation (5.17c) is 

wetted perimeter of channel replaced by Rc D - \'rD/v c . 

=4(1 - p)/oi A ( 5 j 7b ) 5.4. Pressure drop 

In lieu of equations (4.1). (5.15) and (5 17b) NC ^YTT ^ '^T'^ by lhe Cr ° SS - fl ° W 
can be recast into the more familiar form ZnnecSn" ° f .^ Mre| y*™^ Parallel yet 
unconnectmg cylinders may be expressed as [9 16] 

ro-ep^w (5,70 (s,8) 

where C= C(4 P Bi-«0 W tanh(Bi"^. Often, the where is the number of cylinder rows in 
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Fig. 8. Normalized foam performance indices, (a) l u (b) / 2 , as functions of foam density for selected 

values of n. 



f = 0.044 + 



0X m(a/d) l Rp -o.i5 f51 m 



the direction of flow and / is the friction factor. 
(The number of cylinders in the direction perpen- 
dicular to the flow is large so that Ap may be 
assumed uniform in the transverse direction.) For Here, Re max = w//v f and v max = v f /( 1 - dja) is the 
gases,/takes the following empirical relation [16,23] maximum fluid velocity inside the bank. For gases 
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Fig. 9. Prototypical design of a compact multi-chip mod- 
ule cooled by forced convection through open-celled foam. 

flowing across a foam with open cubic cells, the 
pressure differential produced by cylinders of type 
Kand Z may be approximated using equation (5.18) 
and equation (5.19) with N L = 2Lla. The contri- 
bution to the pressure drop from type X cylinders 
and the plate surfaces is more complicated but the 
overall response of such an open-cell foam may be 
characterized by a frictional coefficient, which is 
to be calibrated from experiment. In addition to the 
bulk frictional drag induced by the foam skeleton 
equation (5.18), the fluid flowing past a foam-filled 
duct also experiences a boundary frictional drag at 
the plate wall and a excess pressure loss at the 
entrance region. However, Kaviany [9] showed that 
the total pressure drop experienced by the flow 
across the foam-filled duct is dominated by the bulk 
frictional drag, hence 



(5.20) 



As a reference, the pressure drop of fully-deveioped 
laminar flows in a foam-free duct having length L 
and smooth surfaces is [16] 

&P = (l2L/b)p f v}/Re b (5.21) 

where Re b = V{ b/v r and v ( is the mean flow velocity 
Equations (5.19) and (5.20) suggest that 
Ap^^a 3 " m with m = 1.85 for foams having 
cubic unit cells. Experimental data on metal foams 
is scarce, but for open-celled ceramic foams it is 
found that m has a weak dependence on cell size 
with m = 1.82, 1.88, 1.91 and 1.93 corresponding 
to foams of size 50, 30, 20 and 10 ppi 
respectively [24]. For ERG foams, a smaller expo- 
nent (m - 1.60 both for 20 and 30 ppi foams) has 
been measured [13]. In general, the experimental 
data may be expressed using the quadratic 
interpolation [7J: 

to obtain the Forchheimer coefficient c F and the 
permeability K. For creeping flows with low vel- 
ocities, equation (5.22) reduces to Darcy's law. 

5.5. Optimum foam design 

Ideally, foams for many heat exchanger appli- 
cations would simultaneously maximize the amount 
of heat transfer and minimize the pump power 
needed to force the fluid through. Using 
equation (4.2) and noticing that tanh(2Bi l/2 6/^ % 1 
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in most applications, we can rewrite the overall 
heat transfer coefficient h of equation (5.15) in the 
form 



n/2 



nU2 



which, from equations (5.15) and (5.24), becomes 

(\-dlaf 
/l=P (a/d-\f Ml5 



(a/d--\) n/2 
(5.23) 



where C«0.3 and m%0.6 if a/cf>3. Similarly, we 
can regroup the pressure drop ^p of equation (5.20) 
into 

[ 0M(a/d) 
x 0.044 + 



(a/d-\) 



0.43+1.1 3 (rf/a) 



(5.24) 



where it is assumed that x is independent of ajd. 
The predicted overall heat transfer coefficient, h 
and pressure drop per unit length, Ap/L, are plotted 
on Fig. 7 as functions of relative foam density p 
and fluid velocity v r for air at 20°C flowing across a 
duct filled with 20 ppi (a - 1 mm) ERG aluminum 
foam. These results are generally in agreement with 
the measurements by Schlegel et al. [24] on a cer- 
amic foam except that its overall heat transfer coef- 
ficient is two orders of magnitude smaller than that 
of an ERG metal foam. 

Whilst a foam having the smallest cell size and 
highest relative density maximizes heat transfer, it 
unfortunately also maximizes pressure drop 
(Fig. 7(b)). The predicted minimum pressure drop 
occurs at a relative density roughly equal to 0.09 
for ERG foams, which is close to that measured by 
ERG. The heat transfer performance of the foam 
may be best characterized by a non-dimensional 
index /, combining both Ji and Ap, namely, 



/i = 



QAp 



(5.25a) 
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Prototypical design of a multi-layered counter-cur- 
rent heat exchanger. 



where C x is a parameter (with dimensions) indepen- 
dent of relative foam density, given by 

/ V3^ \ / y /2+0 -' 5 (A,/£)(CPr 1 '*;. f /;,«) l/2 

(5.26) 

For selected values of n Fig. 8(a) plots the index of 
foam performance, hC u as a function of relative 
foam density p. Clearly, there exists an optimum 
foam density, p 0 pt,mum, which would maximize f x C x . 
The value of p optim * m slightly increases as n is 
increased, but generally falls in the range 0.35 - 0.4. 
The cell size, a, does not affect p op timum as well as 
the behavior trend shown in Fig. 8 but, at a given 
density, foams with larger cell sizes perforin better 
than those having smaller cells in size (/i-£f" /2 + 0 1 )- 
The experimental measurements of Schlegel et 
al. [24] on open-celled ceramic foams appear to sup- 
port this latter conclusion. The optimal spacing of 
parallel plates or cylinders in forced convection is 
studied by Bejan [25, 26]. 

If the weight of an heat exchanger is cause for 
concern in applications such as airborne equip- 
ments, the appropriate index for performance scal- 
ing may be taken as h = hlP< Here, the best 
performance is achieved by a foam with p — 0 hav- 
ing the largest cell size a (Fig. 8(b)). 

6. APPLICATIONS 

The model presented in previous sections can be 
used to analyze a variety of heat exchanger appli- 
cations where the large surface density of open- 
celled metal foams is required. Two heat exchanger 
designs are studied below, one being targeted as the 
heat sink for next generation multi-chip modules 
and the other as the air-fuel exchanging system on 
civil/military aircrafts. 

6.1. Heat sinks for power electronic devices 

Heat dissipation in high power electronics poses 
serious challenges for the integration of materials 
selection with thermal design, circuit design and 
manufacturing technology [6,27,28]. The objective 
is to select materials that enable the Si chip to oper- 
ate with high power density (up to 10 7 W/m 2 ), while 
maintaining its temperature below that needed to 
ensure acceptable reliability: usually 90°C 

As an illustration, consider the prototypical de- 
sign of a compact multi-chip module (MCM) envi- 
sioned in Fig. 9. The chips are placed on two plates 
located at z = 0 and z = b that sandwich an open- 
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celled foam, with cooling achieved by forced con- 
vection through the foam. Both plates are taken to 
be good dielectrics and conduct heat well such that 
they are essentially at the same temperature as that 
of the chips, 7V The module is thermally insulated 
at the top and bottom by protective covers (not 
shown in Fig. 9) and, in addition, no heat is lost in 
the ^-direction. A similar design of package has 
been analyzed by Lu et al [6] wherein the foam is 
replaced by a continuum substrate and cooling is 
achieved around the edges x = 0 and x = L. It is 
found that, compared to other package designs con- 
sidered in Lu et al. [7], this particular design suffers 
from inferior heat dissipation, even when high ther- 
mal conductivity substrates with high heat transfer 
strategies, such as water cooling, are used 
However, as elucidated below, the compact design 
of MCM seems more likely to be successful if a 
foam with good heat removing capability is used as 
the substrate, coupled with air or water cooling. 

With a denoting the distance between the centers 
of two neighboring chips (Fig. 9), the amount of 
heat dissipation per chip in a module of length 
L = / is 
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excellent substrate material for the next-generation 
chips: indeed, A1N is currently being used by 
Rockwell International to develop novel high tem- 
perature, wide bandgap solid state power 
devices [30]. Finally, but not least importantly the 
cooling strategy of Fig. 9 dictates that "heat spread- 
ing" in the substrate emphasized in Lu et al. [7] is 
unnecessary, enabling high chip-level power den- 
sities at large chip densities with a * chip size. 



e c = 2A s Bi'/V'/ 2 Ar(^)tanh(^j (6 .l) 

where AF = F c - T 0 . For both water and air cool- 
ing, Q c as a function of p is plotted on Fig. 10 for 
/ s = 200 W/(m K), a = 1 mm (20 ppi), I = 2 cm 
Ar- 70°C, and b/d - 100. It is seen that, even 
with forced air convection in a foam with relative 
density of 0.2, heat is efficiently dissipated from the 
multi-chip module «2 C = 70W if air cooling and 
490 W if water cooling), indeed making it a very 
attractive design for micro-electronic devices due to 
its compactness. 

The selection of a substrate material for holding 
the chips is constrained by many other factors: its 
thermal expansion must be compatible with the Si 
chip to minimize the thermally induced stresses, its 
dielectric constant must be small, should be light 
weight, and it must not contain any toxic ingredi- 
ent. The materials that satisfy these requirements 
include aluminum nitride (AIN), beryllium, alumi- 
num reinforced with particulate silicon carbide and 
berylhum/silicon carbide composites [6,29]. Foams 
made of AIN, for instance, are excellent heat con- 
ductors with x s up to 300 W/(m K). Since the ther- 
mal expansion coefficient of AIN nearly equal that 
of Si, the residual stresses induced in the chips due 
to thermal mismatch are very small. The dielectric 
constant of an AIN foam is exceptionally low if 
cooling is achieved by forced air, 



£ = £f + - er)p 



(6.2) 



where e s = 8.5 is the dielectric constant of solid AIN 
and e f =* I the dielectric constant for dry air All 
these reasons suggest that AIN and its foams make 



6.2. Multi-layer heat exchangers 

Whilst heat exchangers consisting of a single 
channel with constant wall temperatures such as 
that shown in Fig. 1(a) are common, practically it 
is much more common to design heat exchangers 
with more complicated geometries to suit the need 
of transferring heat from one streaming fluid to 
another. Figure 1 1 presents the design of a multi- 
layer counter-current heat exchanger for aeronauti- 
cal applications where hot air flowing in one direc- 
tion is used to preheat jet fuel flowing in the 
opposite direction: this particular design is at pre- 
sent being developed at United Technologies [31]. 
(It is well-known that a co-current exchanger is less 
attractive since it always requires a greater length 
than the counter-current one to perform the same 
duty.) Each layer consists of an open-celled foam 
core, thickness b and length L, sandwiched between 
two thin sheets: the foam and the sheets are made 
of the same material, usually a high conducting 
metal such as aluminum or nickel. The layer is infi- 
nitely wide in the v-direction with negligible end 
effects on heat transfer. For air let m a , p., v a , r a 
and v a denote its mass flow rate per unit wTdth! 
density, kinematic viscosity, temperature, and mean 
velocity inside each layer, respectively. The sub- 
script f will be used to denote the corresponding 
quantities for the fuel. Subscripts 1 and 2 will refer 
to the sections corresponding to inlet and outlet of 
the air layers (or, equivalently, outlet and inlet of 
the fuel layers). The heat exchanger consists of N 
identical layers (only six layers are shown in 
Fig. 11), with hot air flowing in half of them and 
jet fuel streaming in the rest. The midplane of each 
layer (dashed line on Fig. 11) j s thermally insulated 
due to symmetry. Denote by h 0 the overall heat 
transfer coefficient across each sheet. 

The sheet temperature is neither known nor con- 
stant for the case considered in Fig. 1 1, yet one can 
proceed without knowing it, as demonstrated 
below. The air-side heat transfer coefficient, S a , can 
still be calculated according to equation (5.21) 

h a = {*A/2)y/h a Xid (6.3) 

where Bi l/2 b/d > 2 has been assumed. This relation 
also holds for the fuel-side heat transfer coefficient, 
*r, except that the local heat transfer coefficient h[ 
for air across a bank of cylinders is replaced by h ( 
for fuel (cf. equation (4.2)). In many situations, the 



LU et ai: HEAT TRANSFER IN OPEN-CELL METAL FOAMS 



3634 

fluid-related properties, Ji a , h K and h u can be 
safely taken as independent of either temperature or 
position along the sheet. If the sheets that separate 
air and fuel are thin and the fouling factor (due 
either to corrosion products or accumulation of 
various deposits) is neglected, the overall heat trans- 
fer coefficient of the system is given by 

(6.4) 



/io K h( 



If however the sheets are thick or the fouling factor 
is large, extra terms must be added to equation (6.4) 
to allow for the additional thermal resistance. 
Following the method outlined in Section 4.2, an 
enthalpy balance on the element of the exchanger 
shown in Fig. 11 gives the total heat transfer rate 
per unit width, Q, as 



Q = (NL)h 0 AT m 



(6.5a) 



where AT m is the logarithmic mean temperature 
difference defined by 

(r a2 - Tn)-(Tn-T n ) 



in[(r a2 - r f2 )/(r al - 7h)l 



(6.5b) 



If the velocity of air is high at a high temperature, 
then r al *r a2 and n 0 *h »n which case 
equations (6.5aHW reduce to that corresponding 
to a single duct system having a constant wall tem- 
perature. 

The pressure drop in the air layer Ap a follows 
directly from equation (5.26) as 

A/7 a =(2L/a)/aXaPaV am ax < 6 ' 6 ) 

where the maximum velocity in the layer is given by 
2m a 



' Nbp z (l -did) 



(6.7) 



On the fuel side, the pressure drop is calculated 
from equations (6.6) and (6.7) with each sub- 
script a replaced by f. Since the pressure drop 
varies at LPT m b' m with m*1.84 f we can always 
select the best combination of L, N and b to 
lower the pumping power once it becomes unac- 
ceptable. A high pressure drop not only increases 
the operational cost but often also leads to intol- 
erable noises. One of the advantages of designing 
exchangers with foams is that they still offer 
excellent heat transfer coefficient even if the 
mean flow velocity falls below the critical value 
required to maintain turbulent flow in an other- 
wise foam-free duct. Thus, whilst the layer thick- 
ness b is limited by the foam cell size a, the 
restriction placed on the total number of layers 
is much less stringent. Also, the foam itself may 
act as a noise suppressor due to its sound 
absorption capabilities. 
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APPENDIX A 

Heat transfer due to radiation and convection from a heated 
cylinder 

Consider a thin, isolated cylinder of uniform temperature 
T and diameter d. standing in a transverse flow of uniform 
temperature T x { < T). The amount of heat transfer due to 
convection per unit length of the cylinder is 

qh = (nd)h{T-T f ) (A. I) 

while that due to radiation per unit length is approxi- 
mately 

q T = {nd)h I (T- T f ) (A.2) 

Here, h is the average convectivc heat transfer coefficient 
and h T the effective radiative heat transfer coefficient given 
by 



h t = ae{T 2 + Tj){T-r 7» (A.3) 

where a = 5.67 x 10"* W/(m 2 K 4 ) is the Stefan-Bolzmann 
constant and c is the normal total emissivitv of the sur- 
face. For aluminum, t is in the range 0.04-0.2 [16], and 
the maximum service temperature for aluminum-based 
metal foams is about 500 K [2]. For turbulent flows which 
are expected to occur in a bank of cylinders representative 
of metal foams, the convective heat transfer coefficient h 
calculated from equation (4.2) is typically on the order of 
I0--10 W/(nr K), while that of the radiative coefficient h t 
based on equation (A.3) is 10° - I0 1 W/(m 2 K). 
Consequently, heat transfer due to radiation is negligible if 
metal foams are applied in heat exchangers involvine tur- 
bulent flow of fluid. At temperatures above 1000 Eradi- 
ation is expected to play an important role in electrically 
heated thin wires made of refractory metals [32]. 
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